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^ . Abstract 

(D ■ 
pLn I Questions concerning nuclear phase transitions are addressed within the 

QQ ■ transport theory. The issue of thermal equilibrium in the region of the 

^~~i , liquid-gas phase transition is investigated. Signatures of the transition 

to the quark-gluon plasma are looked for. 

>: 

tJ- '. 1 Introduction 

O; 

CN I Under the idealized conditions of thermal equilibrium, the liquid-gas phase- 

Q^ I transition for nuclear matter is a necessity [|^ . However, is there enough time in 

0\ • the reactions to reach the equilibrium? The numerical QCD lattice calculations 

^ ■ exhibit a phase transition to the quark-gluon plasma (QGP). However, can the 

I . excited dense phase leave behind any clear signals in reactions? 

O I We shall try to address the above questions within the transport model 

^ I based on the Boltzmann equation. 

.£;■ 

K> ! 2 Boltzmann Equation 

H ■ 

Cq \ The Boltzmann equation has the same general form relativistically as nonrel- 

ativistically: 

Here, ep is the single-particle energy, 9ep/9p is the velocity, — 9ep/9r is the 
force, and /C™ and /C°"* are the feeding and removal rates, respectively. The single- 
particle energy and momentum, (ep,p), form a 4- vector. Ensuring covariance 
in a model can be aided by using a scalar potential Us dependent on the scalar 
density ps = ^j^ ^^ J dp^ fx- In our model, we also add to the energies 





Figure 1: The deuteron production process run backwards in time becomes 
a breakup process. 

a noncovariant contribution accounting for the Coulomb repulsion, isospin de- 
pendence, and finite-range effects. The single-particle energies are then 



ep = v^^TkTT^F + f^. 



with 
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1 + (e/2.5)-i ' 
^ = ps/po, and 

U^ = Koul + rf V V + C t3 PT (4) 

where the isospin density is px = J2x ^3 P^ ■ 

Besides nucleons and mesons, we consider in our model the lightest clus- 
ters {A < 3) whose phase-space distribution functions follow similar transport 
equations to those for the other particles. An obvious issue with the clusters 
is that of their production. Beyond constituents, a minimum of one addi- 
tional nucleon must be involved in the process, as a catalyst. Let us take the 
deuteron as an example. In the kinetic limit of low rates, the rate for deuteron 
production in three-nucleon collisions is given by the the matrix element for 
production squared convoluted with the (5-functions for the energy-momentum 
conservation and product of statistical factors: 



dpN dpp dpn dp'j^ \M^pN^Nd\2 
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This apparently just transfers the problem of production to the issue of de- 
termining the matrix element for production. Because of the time-reversal 
invariance, though, the production process run backwards in time becomes 
the breakup process, as illustrated in Fig. |l|, and, correspondingly, the two 
processes share the matrix element squared: 
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Figure 2: Liquid-gas phase-coexistence in the plane of pressure and isospin 
asymmetry at T = 10 MeV, for different parametrizations of interactions in 
nuclear matter, after Ref. 0. 

The squared element for the breakup is proportional to the breakup cross 
section 

dcTNd^Nnp 



\M 



Nd^Nnp 1 2 



OC 
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dpdQ 

and, thus, the data on breakup can be used to describe production. The pro- 
duction of y4 = 3 clusters can be handled in a similar manner [01. 



3 Isospin Asymmetry in Low-Density Nuclear Matter 

Under equilibrium, in the phase-transition region for isospin-asymmetric nu- 
clear, the asymmetry is decreased in the denser liquid phase. This is both 
because the interactions play a greater role in the liquid than in the gas phase, 
and they favor the symmetry, and because of the Pauli principle. The en- 
hanced symmetry leaves a pronounced asymmetry in the gas phase. This was 
recently discussed by Miiller and Serot P], and earlier by Barranco and Buch- 
ler m, and by Glendenning ^, and is illustrated in Fig. 0. The difference in 
asymmetry may be as large as by a factor of 3. An obvious issue is whether 
pronounced asymmetries can be seen in reactions and used as a signature for 
the phase coexistence. 
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Figure 3: Relative isotope ratios in Xe-Sn collisions as a function of the veloc- 
ity along the beam in cm 0. Different lines represent results from different 
isotopic Xe-Sn combinations. 



The problem with the gas phase is in its low density. If anything cor- 
responding to that phase is emitted, it will most likely just escape to the 
vacuum. The seemingly only practical way to trap the phase and reach any 
kind of state of phase coexistence is in the neck region between two slow nu- 
clei, with the matter in the nuclei playing the role of a liquid phase. This has 
a further advantage that the gas phase is then well localized in the velocity 
space. Dempsey et al. ^ studied relative yields of different isotopes emitted 
from xenon-tin reactions at 55 MeV/nucleon and found large enhancement 
in the relative abundance of neutron-rich fragments towards midrapidity, as 
illustrated in Fig. ^. 

To verify whether the observed enhancements might signal effects of the 
phase equilibrium, calculations of the reactions have been carried out [0 within 
the transport model with isospin asymmetry in the interactions and in the 
Pauli principle, and with the inclusion of light clusters, as discussed before. 
Figure ^ shows the contour plots for nucleon, proton, and neutron density, 
projected onto the plane of the ^^^Xe -|- ^^^Sn reaction, at different times, 
and the neutron-to-proton ratio along the axis joining the two nuclei. The 
ratio is shown both including the protons and neutrons bound in clusters and 
excluding them. When the simulation is carried without clusters, the results 
for the ratio are similar to those for all protons and neutrons. It is seen that 



regular BUU 
Neutrons Protons • BUU w/clusters 




Ratio 



P 



-10 10 -10 10 -10 10 5 10 15 20 



Z (fm) 



^\) (fm) 
Distance aiong thie neck 



Figure 4: Contour plots for nucleon, proton, and neutron density projected 
onto the plane of ^^^Xe + ^^^Sn reaction, at different times, and the neutron- 
to-proton ratio along the axis joining the two nuclei 0. 

the ratio for all nucleons just slightly exceeds the N/Z ratio of 1.5 for the whole 
system. When the clusters are excluded, the ratio reaches high values, up to 4. 

The reason for the enhancement in the calculation is the deuteron forma- 
tion which robs the remainder of the neck region of protons enhancing the 
asymmetry for the remainder. One can expect that, in reality, alpha particles 
play an analogous role to deuterons leading to the asymmetries observed ex- 
perimentally. On the other hand, the calculations show that the overall isospin 
asymmetry between the gas and liquid phases, due to the favoring of the sym- 
metry in the liquid and expected in the equilibrium, has no time to develop. 
This brings in the general issue of equilibration within the neck region. 

The isospin ratio may be viewed as a convenient variable to assess the equi- 
libration within the phase transition region that, seemingly is not reached. Cal- 
culationally, it is far easier to determine the isospin asymmetry than e.g. pres- 
sure or temperature. On the other hand, one can argue that the isospin equi- 
librates in a slow diffusive process (as e.g. smell, aided, though, in the air 
by convection) while such quantities as pressure or overall density equilibrate 



rapidly in a mechanical process mediated pressure waves. Thus, isospin could 
not be a reliable measure. 

Let us take L as a distance over which adjustments of a thermodynamic 
quantity are to propagate and t as time that the propagation takes. In a dif- 
fusive process, a particle undergoing A^ collisions covers a typical distance: 

L'^ = NX\ (8) 

where A is the mean-free-path. The collision number is related to the net time 
by t = A^ X/v, where f is a typical particle speed. The characteristic time for 
the isospin adjustements is then: 

f-iso T • W J 

Xv 
The time for the density and pressure adjustements, on the other hand, is 

_ L ^ L 

Cs V 

The ratio of the two times, 

''mech ^ 

would have been large, if A were short. However, in the neck region p ~ Po/3 
and thus A ~ 6 fm. At the same time the distance from the nuclei to the 
center of the neck is L ~ 5 fm and we are in the Knudsen limit! There is no 
difference between the times for isospin and density adjustements, both are of 
the order of 30 fm. Even for a Knudsen gas, the equilibration could take place 
through a contact with a thermostat, i.e. nuclei, but that does not happen in 
the simulation. We find that we are not dealing with a phase transition near 
equilibrium. 

4 Model for the Transition to Quark-Gluon Plasma 

Let us now turn to the other important phase transition in nuclear physics, 
to QGP. The phase transition is approached when hadrons increase in num- 
ber, pushing out from their region the nonperturbative vacuum. This can 
make the hadrons lighter as the condensates responsible for hadronic masses 
are removed. The transition is expected to occur when hadrons completely 
overlap. The lattice calculations for baryonless matter point to a phase tran- 
sition at the temperature of Tc ~ 170 MeV corresponding to the pion density 



p^ ~ 0.5 (^) ~ 2po- If the phase transition is associated with the hadrons 
fining up all space, the transition should be expected also in the relatively cold 
matter when hadron density is raised to comparable values due to compression, 
rather than thermal production. The idea is further illustrated with a simple 
quasiparticle model ^j. 

The model 0, based on the hadronic degrees of freedom, produces a de- 
crease of masses in the vicinity of the phase transition, a large drop in masses 
across the transition, and a dramatic increase in the number of the degrees of 
freedom in the transition. At the same time, the properties of ground-state 
nuclear matter are correctly described. The model is specified by giving the 
energy density as a function of particle phase distributions in the form of ki- 
netic part and interaction corrections depending on two densities, scalar and 
vector: 

e = J2jdp4 /""(P) + ^siPs) + e.(p) , (12) 
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where 

and p = Ylx ^^ I ^P f^ip)- '^^^ scalar density counts all hadrons while 
the vector density is just the baryon density. As hadron density increases, 
the hadron masses drop, from (|T^) 



5e 



4 = jjxt::^ = \/p' + (SiPs) m^f + A^ U{p) , (14) 



by a common (for simplicity) factor 5* in our model. In a thermally excited 
system the lowering of the masses leads to the production of more hadrons. 
Instability can occur leading to a phase transition across which the masses 
drop and the number of the degrees of freedom rapidly increases. We adopt 
simple power parametrizations of the mass modification factor and of the vector 
potential to get consistency with lattice calculations and ground-state nuclear- 
matter properties, 

S={1- 0.54 (fmVGeV) p^f , (15) 

and 

l + 6p/po + c(p/po)5/3 ^ ' 

where a = 146.32 MeV, b = 0.4733, c = 51.48. With the transition in the 
baryonless limit at Tc = 170 MeV, we find a corresponding transition at T = 0, 
taking the matter from the baryon density of p ~ 3.5 po to 7po, cf. Fig. ^. 
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Figure 5: Chemical potential and energy per baryon at T = as a function of 
baryon density. The high-density phase transition is indicated with arrows. 

5 Elliptic-Flow Excitation Function 

A possible way of detecting the high-energy nuclear phase transition is through 
the associated changes in the speed of sound. Thus, e.g. in the phase-transition 
region in Fig. ^, the speed of sound in the long wavelength limit (c^ = ydp/de) 
vanishes. 

A sensitive measure of the speed of sound or pressure compared to the en- 
ergy density early on in the reactions is the elliptic flow. The elliptic flow is the 
anisotropy of transverse emission at midrapidity. At AGS energies, the elliptic 
flow results from a strong competition [Q between squeeze-out and in-plane 
flow. In the early stages of the collision, the spectator nucleons block the path 
of participant hadrons emitted toward the reaction plane; therefore the nu- 
clear matter is initially squeezed out preferentially orthogonal to the reaction 
plane. In the later stages of the reaction, the geometry of the participant re- 
gion (i.e. a larger surface area exposed in the direction of the reaction plane) 
favors in-plane preferential emission. 

The squeeze-out contribution to the elliptic flow and the resulting net di- 
rection of the flow depend on two factors: (i) the pressure built up in the com- 
pression stage compared to the energy density, and (ii) the passage time for 
removal of the spectator shadowing. In the hydrodynamic limit, the character- 
istic time for the development of expansion perpendicular to the reaction plane 
is ~ -R/cs, where R is the nuclear radius. The passage time is ~ 2R/{jqVo), 
where Vq is the cm. spectator velocity. The squeeze-out contribution should 
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Figure 6: Calculated elliptic flow excitation functions for Au + Au. The open 
diamonds represent results obtained with a stiff EOS. The open circles repre- 
sent results obtained with a stiff EOS and with a phase transition. 



then reflect the ratio 



70^0 



(17) 



According to ([T7| ) the squeeze-out contribution should drop with the in- 
crease in energy, because of the rise in vq and then in 70. A stiffer equation of 
state (EOS) should yield a higher squeeze-out contribution. A rapid change in 
the stiffness with baryon density and/or excitation energy should be reflected 
in a rapid change in the elliptic flow excitation function. A convenient measure 
of the elliptic flow is the Fourier coefficient (cos 20) = f2, where is the az- 
imuthal angle of a baryon at midrapidity, relative to the reaction plane. When 
squeeze-out dominates, the Fourier coefficient is negative. 

To verify whether the expectations regarding the elliptic-flow excitation 
function are realistic, we have carried out Au + Au reaction simulations [p!0| , 
in the energy range of (0.5-11) GeV/nucleon. The excitation functions cal- 
culated using a stiff EOS with a phase transition (open circles) and a stiff 
EOS with a smooth density dependence are compared in Fig. |^. For low 
beam energies (< 1 AGeV), the elliptic flow excitation functions are essen- 
tially identical because the two EOS are either identical or not very different 
at the densities and temperatures that are reached. For 2 < Eseam ^ 9 AGeV 



the excitation function shows larger in-plane elliptic flow from the calculation 
which includes the phase transition, indicating that a softening of the EOS 
has occured for this beam energy range. This deviation is in direct contrast 
to the esentially logarithmic beam energy dependence obtained (for the same 
energy range) from the calculations which assume a stiff EOS without the 
phase transition. Present data on elliptic flow from EOS, E895, and E877 
Collaborations [|ll|] point to a variation in the stiffness of EOS in the region of 
~(2-3) GeV/nucleon, corresponding to baryon densities of ~ 4po- 

6 Conclusions 

To summarize, the observed large enhancements in the yields of neutron-rich 
clusters in the neck region of reactions seem to be due to the deuteron and 
alpha production. Phase transition to QGP at high T in baryonless matter 
implies a transitional behavior at T = and high baryon density. Elliptic flow 
measurements can decide about the presence or absence of the QGP transition 
with the rising density. 
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